Experimental Investigation of Resonant Activation 
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We experimentally investigate the escape from a metastable state over a fluctuating barrier of a 
physical system. The system is switching between two states under electronic control of a dichoto- 
mous noise. We measure the escape time and its probability density function as a function of the 
correlation rate of the dichotomous noise in a frequency interval spanning more than 4 frequency 
decades. We observe resonant activation, namely a minimum of the average escape time as a func- 
tion of the correlation rate. We detect two regimes in the study of the shape of the escape time 
probability distribution: (i) a regime of exponential and (ii) a regime of non-exponential probability 
distribution. 



The thermal activated escape of a particle over a po- 
tential barrier in a metastable state is a classical prob- 
lem since the seminal work of Kramers Q . It occurs in a 
wide variety of physical, biological and chemical systems 
When the height of the potential barrier is itself 
a random variable the dynamics of the escape from the 
metastable state is affected by the statistical properties 
of the random variable controlling the barrier height. In 
the presence of barrier fluctuations rather counterintu- 
itive phenomena may occur. A classical example is the 
phenomenon of resonant activation ||. The signature 
of the resonant activation phenomenon is observed when 
the average escape time of a particle from the metastable 
state exhibits a minimum as a function of the parameters 
of the barrier fluctuations. The original work of Doer- 
ing and Gadoua triggered a large amount of theoretical 
activity Pd^- Analogic simulations of resonant activa- 
tion have also been performed in a bistable system in the 
pre sence of multiplicative Gaussian or dichotomous noise 

This letter aims to answer the following questions - 
(i) is resonant activation experimentally observable in a 
physical system? (ii) what about the probability density 
function of escape time? We answer both questions by 
investigating the escape from a metastable state over a 
fluctuating barrier in a physical system. The physical 
system is a tunnel diode biased in a strongly asymmetric 
bistable state in the presence of two independent sources 
of electronic noise. Specifically we have a dichotomous 
noise source controlling the metastable potential and a 
Gaussian noise source mimicking the role of a thermal 
source. In this physical system we observe resonant acti- 
vation and we measure and characterize the probability 
density function (pdf) of escape time by finding two dis- 
tinct regimes. 

Our physical system is the series of a resistive network 
R with a tunnel diode in parallel to a capacitor (in our 
case the sum of diode capacitance and input capacitor of 
the measuring instrument). This system is rather versa- 



tile and allows investigate the dynamics of a bistable or 
metastable system in the presence of noise. It has been 
used to investigate stochastic resonance |l8| and noise 
enhanced stability ]19| ], 

The equation of motion of such system is 
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where R is the biasing resistor of the network, C is the 
capacitance in parallel to the diode, V rrns is the ampli- 
tude of the Gaussian noise v n (t), Vb is the biasing voltage 
and I(v) is the nonlinear current-voltage characteristic of 
the tunnel diode. The effective associated potential is a 
bistable potential. We control the parameters of our elec- 
tronic network (Vb and R) in such a way that one of the 
two wells is much deeper than the other. By setting such 
strong asymmetry between the two wells we essentially 
deal with a metastable state. In fact, the probability that 
the system once escaped goes back into the metastable 
state during the experimental time is negligible. By using 
a digital electronic switch we vary the value of the scries 
resistance of our metastable system between the two val- 
ues R+ and R-. When the system switches from R + to 
R- two effects take place. One concerns the variation of 
the height of the barrier of the metastable state whereas 
the other consists in a slight change of the value of intrin- 
sic time scale of the system (the RC term of Eqs. (1) and 
(2)). In our experiment we use a germanium tunnel diode 
1N3149A, which has nominal peak current I p = 10.0 rnA, 
peak voltage V p — 60 mV, valley current I v = 1.3 mA, 
valley voltage V v = 350 mV and a rather short switching 
time (less than 50 nsec). We perform our experiments by 
choosing V b = 8.95 V, R+ = 1100 fl, i?_ = 1080 and 
C = 45 pF. The noise v n (t) is obtained starting from 
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a digital pseudo-random generator. The noise voltage is 
a stochastic Gaussian noise synthesized by a commercial 
source (the noise generator DS345 of Stanford Instru- 
ments). The Gaussian noise v n (t) is added to the bias 
signal Vb by an electronic adder based on a low-noise 
wide band operational amplifier. At the output of the 
operational amplifier we measure the noise v n (t). It is a 
Gaussian noise characterized by a spectral density which 
is flat at low frequency (/ < 1 MHz), having a moderate 
(approximately 7 dB) increase in the region (1 < / < 3 
MHz) and quickly decreasing after the cut-off frequency 
(/cni-o//=4.6 MHz). By defining the correlation time of 
the Gaussian noise as the time at which the normalized 
autocorrelation function assumes the value 1/e, we mea- 
sure r n = 68 ns. In addition to the Gaussian noise source, 
a dichotomous noise source is also present. The dichoto- 
mous noise source is obtained by using the commercial 
chip MM5437 of National semiconductor. An external 
clock drives the dichotomous noise source and allows us 
to control the noise correlation time over a wide range. 
The normalized autocorrelation function of this pseudo- 
random noise is linearly decreasing from 1 to zero in one 
clock period l// c . After this time the autocorrelation 
function is equal to zero within the experimental errors. 
By defining the correlation time as before, we have that 
7 = T(j l ~ 1.3 x f c . In our experiments the correla- 
tion rate 7 is varied in the range from 7 m i„ = 18 Hz to 
Imax — 316228 Hz. This is an experimental interval con- 
sisting of more than 4 frequency decades. To investigate 
a so wide range of frequency we have to overcome two ex- 
perimental conflicting constraints: (i) we are forced to set 
the time constant of our system r s = RC to a low value 
satisfying the inequality ^max *C 1/t s and (ii) we need 
to use a high value of t s to maintain the ratio t u /t s as 
low as possible to conduct our experiments in the "white 
noise" limit of v n (t). The best compromise we find is to 
set t s < R + C = 49.5 ns. With this choice l/r s rs 507 maa; 
and t„/t s rs 1.40. In other words we guarantee the in- 
vestigation of the resonant activation phenomenon in a 
wide range of frequency by performing our experiments 
in a regime of moderately colored noise. 

Under computer control, we measure, for each value of 
the correlation time of the binary noise, 5 x 10 3 escape 
times from the metastable fluctuating state. In each sta- 
tistical realization, the system is put into the metastable 
state by an electronic switch. We define the escape time 
T as the time interval measured between the setting of 
the system in the metastable state (t = 0) and the cross- 
ing of Vd of a voltage threshold. The selected voltage 
threshold (vd — 0.4 V) ensures that the system is quite 
far from the starting well. The exact value of the thresh- 
old is not a determinant parameter because the escape 
from the well is fast (of the order of t s ) . From the set of 
measured escape times T, we determine the average value 
< T > and the distribution P{T). The first measurement 
is devoted to determine the Kramers rates as a function 
of the noise amplitude V rms for the two metastable states 
in the absence of the dichotomous noise modulation. In 



our experiments we vary V rms in the interval from 0.816 
V to 1.00 V and we observe that the logarithm of the 
average escape time is a linear function of the inverse of 
the noise intensity 1/V? ma for the two states "+" and 
"-" of the metastable state. In other words, the mea- 
sured values lay in a straight line verifying the Kramers 
law in both cases. The two metastable states are charac- 
terized by different slopes of the two lines fitting exper- 
imental data. These experimental observations confirm 
that the shape of the metastable potential is different 
in the two states "+" and "-" . The observation of a 
Kramers law is in agreement with the theoretical pre- 
diction of the escape from a metastable state both in the 
absence and in the presence of colored noise. Specifically, 
in a regime of colored noise v n (tY the expected functional 
form is still Kramers-like J20|]21|] and can be written as 



< T >= C exp [AU m T^ 



221 . Here C is a prefac- 



tor and AU m is the measured potential barrier height as- 
sociated to a system characterized by a time constant r s 
and evolving in the presence of a colored noise with cor- 
relation time r n . Theoretical studies predict that AU m 
is a function of r„/r s 20-22|. We verify that we are in 



a regime of colored noise by performing an experimental 
test 

We perform our investigation of the average escape 
time as a function of the correlation rate of the dichoto- 
mous noise for two different values of the amplitude of 
the Gaussian noise. The chosen values are V rms = 0.816 
V and V rrns = 0.852 V. In the absence of switching be- 
tween the two states, the average escape times from state 
"+" and "-" are < T+ >= 0.030 s and < T_ >= 0.0021 
s when V Tma = 0.816 V and < T+ >= 0.0082 s and 

< T_ >= 0.00081 s when V rms = 0.852 V. In the pres- 
ence of switching between the two states the behavior 
is different. In Fig. 1 we show the average escape time 
as a function of 7 in a log- log plot. We also indicate 
the value of the Kramers rates /i+ = 1/ < T + > and 
/i_ = 1/ < T_ > for the two sets of experiments. In 
both cases the distinctive characteristics of resonant ac- 
tivation are observed ||] , the average escape time initially 
decreases, reaches a minimum value and then again in- 
creases as a function of 7. Specifically for values of the 
correlation rate 7 less than /i + the average escape time 
approaches the value < T >= (< T_ > + < T + >)/2 
predicted by resonant activation theory ^J6],^3|. We ad- 
dress this regime as the regime of slow dynamics of the 
barrier height. For values of the correlation rate approx- 
imately ten times higher than the minimum value 
of the average escape time is observed. Here resonant 
activation is fully effective. In this regime, the value of 
the average escape time is in both cases approximately 
10% less than the one expected from resonant activa- 
tion theories in the kinetic approximation PJl3]] , namely 

< T >= ((/z+ + /i_)/2) _1 . We have not a definitive 
explanation for this discrepancy. One possibility is that 
this effect reflects the fact that in our experimental set- 
up our system switches between two states in a regime 
of colored Gaussian noise rather than in the regime of 
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"white noise". A second one is that this behavior is re- 
lated to the fact that our system switches between two 
states which in addition to a different potential barrier 
height are also characterized by different shape of the po- 
tential and different time constant. For values of 7 higher 
than f kHz the average escape time starts to increase as 
a function of 7 and approaches the value associated with 
an effective potential characterized by an average barrier 
at the highest value of 7. This is the regime of fast dy- 
namics of the barrier height. A comparison of the two 
sets of experiments of Fig. f also shows that the region of 
resonant activation becomes more flat when the noise in- 
tensity decreases. For our experiments the barrier heights 
are kept constant in the two series of experiments. This 
observation is complementary to the theoretical observa- 
tion that the resonant activation region becomes more 
flat when the barrier height is increased while the noise 
intensity is kept constant Jl^ , p^] . 

For each value of the noise intensity V^ ms and for each 
value of the correlation rate we measure P(T), the pdf of 
the escape time. In Fig. 2 we show all the pdfs measured 
with V rms = 0.816 V in a 3-dimensional semi-logarithmic 
plot. The regimes of slow dynamics, resonant activation 
and fast dynamics of the dichotomous noise are clearly 
observed moving from the left to the right. Concerning 
the shape of the pdf, our results show that two distinct 
regimes are present. In particular for values of 7 greater 
than P(T) is well described by an exponential decay 
function P(T) = exp [-77 <T >]/ <T > whereas in 
the opposite regime the pdf is non-exponential. Hence, in 
the resonant activation regime, the pdf has an exponen- 
tial shape. This experimental observation is in agreement 
with the numerical observation that the system escapes 
preferentially through the state with the lowest barrier 
at the resonant activation In other words, at the res- 
onant activation the system approximately experiences a 
single potential barrier and this ends up into an exponen- 
tial pdf. Exponential pdfs are also observed for higher 
values of 7 in the fast dynamics regime. The experimen- 
tal detection of an exponential pdf for high values of 7 is 
in agreement with the simple description that the system 
experiences an effective potential with an average barrier 
for the highest values of the correlation rate. 

The final investigation of our study concerns the de- 
gree of non-exponential behavior observed when 7 < 
To quantify the difference between the real pdf and an 
exponential pdf we measure the standard deviation of 
the average escape time. For an exponential pdf the two 
observables coincide, hence a difference in these observ- 
ables manifests a non-exponential behavior of the pdf. In 
Fig. 3 we show the average escape time together with its 
standard deviation as a function of 7 measured by set- 
ting V rms = 0.816 V . The two regimes are clearly seen. 
For high values of 7 the pdf is well described by an ex- 
ponential function, here the average escape time and its 
standard deviation essentially coincide. When the corre- 
lation rate is less than /i_ a deviation from the exponen- 
tial form starts to emerge and becomes more pronounced 



for lowest values of 7. In particular, for very low values 
of the correlation rate (7 < /i+ < /Lt_) the pdf assume 
the form of a double exponential 

P(T) = 2< l r+> cxp [-T/ < T + >] + 

2< y_ > exp[-T/<T_>]. (3) 

The inset of Fig. 3 shows P(T) measured when 7 = 13 
Hz in a semi- logarithmic plot. The shape is clearly non- 
exponential. It is well approximated by Eq. (3) when 
the measured values of < T + > and < T_ > are used 
(solid line in the inset). In this regime the correlation 
rate is smaller than Kramers rates of both states. This 
implies that the system essentially starts and ends each 
realization of the stochastic process in one of the two 
states with probability 1/2 and the overall pdf is given 
by Eq. (3). 

In this letter we present an experimental study de- 
tecting resonant activation and we verify some theoret- 
ical results obtained in model systems for limit regimes 
of the correlation rate. Resonant activation is observed 
in spite of the fact that the Gaussian noise mimicking 
the presence of temperature is colored instead of being 
"white" . In our opinion, our result shows that the res- 
onant activation phenomenon is pretty robust and may 
be observed in a variety of physical systems. By inves- 
tigating the pdf of the average escape time, we detect a 
non-exponential shape of the pdf in the regime of low val- 
ues of 7. This finding suggests that systems characterized 
by a non-exponential escape time probability distribution 
could be simply modeled in terms of metastable system 
with a fluctuating barrier. 
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FIG. 1. Average escape time as a function of the correlation 
rate of the dichotomous noise. The amplitude of the Gaussian 
noise is set to Vrms = 0.816 V (squares) and Vrms = 0.852 

V (circles). For both sets of data, dashed and dot-dashed 
lines indicate the Kramers rate of the two states "+" and "-" 
respectively. Three regimes are clearly seen: (i) the regime 
of slow dynamics between states (7 < pi+); (h) the regime of 
resonant activation (10 3 < 7 < 10 5 Hz) and (iii) the regime 
of fast dynamics between states (7 > 10 s Hz). The solid 
lines indicate the asymptotic behavior expected for slow dy- 
namics whereas the long-dashed lines indicate the asymptotic 
value expected from the kinetic approximation. The regime 
of resonant activation is broader for the value V rms = 0.816 

V (squares). 




FIG. 2. Overall summary of the probability density func- 
tions of escape time measured by setting Vrms = 0.816 V. The 
z axis is logarithmic. From left to right the regimes of slow 
dynamics, resonant activation and fast dynamics are clearly 
seen. In the regimes of resonant activation and fast dynam- 
ics the probability density function is well described by an 
exponential function whereas in the slow dynamics regime a 
non-exponential shape emerges. 
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FIG. 3. Average escape time (black circles) and its stan- 
dard deviation (white circles) as a function of the correlation 
rate of the dichotomous noise. The amplitude of the Gaussian 
noise is set to V rm s = 0.816 V. The difference between the two 
observables quantifies the degree of non-exponential behavior 
of the pdf. Two regimes are observed: (i) the regime of val- 
ues of 7 > fi- where an exponential pdf is observed within 
experimental errors and (ii) the opposite regime (7 < /i_) 
where a non-exponential behavior is observed. The degree of 
non-exponential behavior is maximal in the regime of slow dy- 
namics. In the inset we show the probability density function 
(circles) of the escape time obtained by setting 7 = 13 Hz. 
The pdf is well approximated by Eq. (3) with < T + > and 
< T_ > measured in the absence of switching. This function 
is shown as a solid line in the inset. 
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